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Abstract 

We study the stability of a special form of D brane embedding which is regarded as 
a bound state of D n and anti-D n -brane embedded in a 10D supergravity background 
which is dual to a confining gauge theory. For D5 branes with U(l) flux, their bound 
state configuration can be regarded as the baryonium vertex. For D branes of n = 6 and 
8 without the U(l) flux, their bound states have been used to introduce flavor quarks 
in the dual supersymmetric Yang-Mills theory. In any case, it would be important to 
assure that they are free from tachyon instability. For all these cases, we could show 
their stability with respect to this point. 



gour okuOdont aku . f it . ac . j p 
2 nakamura@sci . kagoshima-u. ac . jp 
3 f toyoda@f uk .kindai .ac.jp 



1 Introduction 



It would be a challenging problem to make clear the deep structure of mesons, baryons 
and other hadrons which are the bound states of quarks. The quantum chromody- 
namics is well defined to solve these problems, but it is very difficult to see the non- 
perturbative properties of this theory. A very promising approach in this direction has 
been performed in the recent holographic approaches. According to the quark model, 
the baryons are composed of N c quarks combined at a vertex point, which would 
however not a point but has a structure and the vertex with a structure has been 
expressed in type IIB theory by D5-branes wrapped on S 5 in AdSs x S* 5 space-time 

miSlSHllSISISUfflmHSHSIIl based on the string/gauge theory correspondence 
[T5| IT6l ITT] . And the baryon has been constructed as a system of fundamental strings 
(F-strings) corresponding to the quarks and this vertex jT8], [191 EH EH E2J E2J EI]- 

In such models, the F-strings are dissolved as a U(l) flux in the D5 brane, and their 
remaining parts flow out from S 5 through one (or two) cusp(s) on the surface of the 
D5 brane to the outside of the extended vertex. The outside fluxes are separated as N c 
free strings which are connected to the same number of separated quarks. The vertex 
has been embedded on S 5 in 10D bulk, and its shape could be deformed due to the 
dissolved U(l) flux and also by the 10D bulk backbacground which is dual to the YM 
gauge theory in the confinement phase. This deformation can also be seen in our real 
3D space as a string like object [31 IU EJ [H] . 

As a special deformed solution of the embedded D5 brane, we have found a config- 
uration which corresponds to a bound state of D5 and anti-D5 brane [32]. And this 
configuration has been regarded as a baryonium (a bound state of baryon and anti- 
baryon) vertex since this has two cusps where a definite number (N q ) of U(l) fluxes 
enter from one side and the same number of fluxes come out from the other cusp. So 
this vertex connects N q quarks and the same number of anti-quarks. So this state can 
be considered as a meson with a special vertex, and it is called as baryonium. So our 
holographic approach have a chance to resolve the long standing problem whether this 
kind of state could be realized or forbidden. So it is an important issue to study the 
stability of our baryonium configuration. 

Previously, we have resolved this problem from the view point of the energy, and we 
find that there is a minimum point of the classical action for an appropriate solutions of 
the baryonium. However, in general, we should be aware of an instability coming from 
tachyons which are anticipated in this kind of configurations, the bound state composed 
of D and anti-D branes. This time, we examine this point by studying the spectrum of 
the fluctuation for the field living on the D5 brane to search for the tachyonic modes. 

We also study the stability of the other cases of the D and anti-D brane bound state 
used in the other holographic gauge theory. 

In the next section, we give our model and D5-brane action with non-trivial U(l) 
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gauge field, and the baryonium solutions are briefly reviewed. In section 3, we give 
a method to study the eigenvalue of the fluctuation modes and the stability of the 
baryonium is discussed. The stability of the bound states of D8 and anti-D8, and D6 
and anti-D6 branes, which are considered to introduce fundamental quarks, are also 
studied by the same method in section 4. Finally in the final section, we summarize 
our results and discuss related problems. 



2 D5/D5 and Baryonium 

The D5-D5 configurations have been proposed as baryonium vertex. The configura- 
tions have been given as solutions of the equations of motion of the D5-brane which 
is embedded as a probe in a supersymmetric lOd background of type IIB theory. The 
dual theory of this background corresponds to the confining gauge theory, then the 
configuration given here can be considered as a bound state of baryon and anti-baryon. 
How to obtain this solution is briefly reviewed in the following. 



2.1 Bulk background and D5 brane action 

We consider the following supergravity background solution 

ds\ = e*' 2 (^Ti^dafda? + ^dr 2 + R 2 dnl)j , (1) 

which is written in string frame. And the dilaton $ and the axion x are given as 

e* = l + ^, X = -e-* + Xo, (2) 
with self-dual Ramond-Ramond field strength 

(7( 5 ) = dC {A) = 4i? 4 (vo\{S 5 )d6 l A...Ad6 5 - ^dt A . . . A dx 3 A dr^j , (3) 

where vo^S* 5 ) = sin 4 ^ivo^S* 4 ) = sin 4 6\ sin 3 62 sin 2 #3 sin 64. 

This solution, ([I])-©, is useful since the confinement of quarks are realized due to 
the gauge condensate q = (F 2 U ) [2S1I2I]' which is given by the coefficient of 1/r 4 for the 
asymptotic expansion of the dilaton at large r. And furthermore, M=2 supersymmetry 
is preserved in spite of the non-trivial dilaton is introduced. We can assure through 
the Wilson loop that q 1 ^ 2 is proportional to the tention of the linear rising potential 
between the quark and anti-quark [27]. In the present case, q is essential to fix the size 
of the baryonium and stabilize it energetically. We notice that the axion \ corresponds 
to the source of D(-l) brane and it is Wick rotated in the supergravity action. This is 
necessary to preserve the supersymmetry. 
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In the next, we introduce the probe D5 brane, and its action must include N c dis- 
solved U(l) fluxes in it. The D5-brane action is thus written as by the Dirac-Born-Infeld 
(DBI) plus WZW term H] 

S D5 = -T B J d^e-^-det (g ab + F ab ) + T 5 J d 6 £A (1) A C {5) , (4) 

gab = d a X fJ, d b X u G fJ-l/ , C ai __, a5 = <9 ai X w . . . d aB X B C )Jll ... l j l5 . 

where T 5 = 1/ (g s (2%) 5 l s 6 ) is the brane tension, and U(l) worldvolume field strength is 
expressed by F ab = dA^ with F ab /2iia' = F ab = dA^y And denotes the induced 
five form field strength. . 



2.2 Equations of motion and Baryonium solution 

The D5 brane is embedded in the world volume £ a = (t, 9, 9 2 , . . . , 6*5), where (6 I 2 , . . . , 6* 5 ) 
are the S" 4 part with the volume of SI4 = 87r 2 /3, where we set as B\ = 6. Then, we 
restrict our attention to 5*0(5) symmetric configurations of the form r(8), x(9), and 
A t (6) (with all other fields set to zero). In this case, the above action is written as [8] 

S = T 5 fi 4J R 4 J dt d6{- sin 4 9 \fe* (r 2 + r' 2 + (r/R) A x 12 ) - F t 2 e - F te D}, (5) 

where the WZW term is rewritten by partial integration with respect to 6, and ^4 = 
87r 2 /3 is the volume of the unit four-sphere. The factor D{9) is defined by 

d D = -4sin 4 fl , (6) 

and is related to F t g by the equation of motion for A t given below. We call this D as 
displacement. It is given by solving (jHJ) as follows , 



D(u,6) = 



3 3 

- (utt — 6) H — sin 6 cos 9 + sin 3 9 cos 9 

— — 



(7) 



Here v is an integration constant, defined in the range of < v < 1. The meaning of v 
is described in details in [U |8], and we can briefly express it as follows. The D5 brane 
as a baryon vertex has two cusps where U(l) fluxes come out the D5 form. The total 
number of the flux is N c , and they are separated to N c (l — v) and N c v to each cusp 
point. Namely, the meaning of v is then the ratio of this flux separation. 

Equations of motion 

Here we comment on the equations of motion which we used actually in obtaining 
the baryonium state and the formal one which can be obtained from the linear terms 
of the fluctuations of the fields on the brane. For the latter case, starting from (j4j) or 
the equations of motion for r(9), x(9) and A t (9) are obtained as 



sin 



4 ^( e */ 2 ^ ( o))-a e [sin 4 ^4^r'| = 



*/2 
(0) 



3 



r 



4 e $/2 N 



djsin 4 # ? =x I = 



p -*/2 - 



*(0) 



where r' = dgr(9), x' = dgx{9) and 



^(0)=(r 2 + (r') 2 + ^(x') 2 -e-^j . (9) 
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The latter two equations are solved as 



sin 4 9 x = h 



r 



4 g $/2 



R\K, 



sin 4 9-=A' t - D = 0, (10) 

V^(0) 

where h is a constant and identified with the one used in our previous paper by the 
same notation [32]. Using these, we obtain the equation for r(9), which will be solved 
numerically due to its complicated form. 

In order to solve in more understandable way, it is convenient to use the equations 
of motion considered before in |32j, and they are given in the Appendix A. It is easy 
to prove the equivalence of the solutions of the above equations (jHJ) and the one given 
in [32]. 

We notice here that the parameter h given in the above equations (|T0|) is equivalent 
with the one of equation (163|) in the Appendix A. Secondly, we briefly review how our 
solution can be identified with the bound state of D5 and anti-D5 branes according to 
the previous work. 

Approximate solution: 

We firstly show how the region of 9 in the baryonium solutions is restriced by using 
h. From Eq. ( 165|) . we find 



then it leads to 



> . (12) 

For prototypical small energy solution, r takes its minimum r ~ r at the mid point 
x = 0. Expressing as 9\ x= o = 9 at this mid point, the value of 9q is given by f fT2|) as 
the solution of the following equation, 

r ^T-V v (9 ) = h* . (13) 




This gives two solutions, 9 and 0\ with 9 > 9±, as shown in the Fig. [TJ Then we find 
that the solutions are restricted to the region < < 0± or 6q < 6 < n , and they are 
identified with two different baryoniums. 

At the mid point, x = 0, 6 = 9q and r = ro, which are related as 

ro = lmH ' (14) 

Near this point, we can assume f ~ and r ~ r , and we obtain from (TTTj) 



9 = + qy/V v (6) - V v %) (15) 

then it is solved as 

x(6) = ± / d0 . (16) 

This solution is symmetric with respect to x = axis in x-6 plane. The important 
point of this approximate solution is that the solution runs from x = to two opposite 
directions, however they are going to the same pole on S 5 but with different values of 
x. These two parts are considered as D5 brane and anti-D5 brane, and they can be 
connected at x = by a continuous flow of the U(l) flux in the D5 brane. The similar 
behavior has been seen in other D brane and anti-D brane bound states j28j [291 130] - 
We regarded this as the baryonium since the fluxes at the two end points of this 
configuration have the same magnitude and opposite directions. 

Numerical solution 
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In the next, in order to see the behavior of this solution far from 6q, we must solve 
the exact form of equations as shown in [32]. An explicit example of the baryonium 
solution for v = 0.3 is shown in the Fig. [2j This solution gives a minimum of U and 
there is no tachyonic mode around this configuration. So this configuration is stable 
as shown below. 



Fig. 2: A 3D representation of D5/D5 solution for R = 1, q = 0.3, v = 0.3 at h = — .67 and 
9q = 1.42, where U takes its minimum and tachyon free. 



3 Stability of the baryonium vertex 

Here we study the stability of the baryonium solution obtained above. Previously, we 
have considered this problem from the energy density U of the baryonium configura- 
tions to find its minimum as a stable solution. And we have found a minimum for 
v = 0.5 as such an example. It has been shown through the contour graph in the h — L 
plane, where L is defined as 



The baryonium solutions are considered as the bound states of D5 brane and anti-D5 
brane. Then they should be separated by a definite distance to evade tachyon modes 
on the D5 branes. The distance between these branes is characterized by L, and the 
U m i n is found in the region of L 2 /R A > — This implies that the tachyonic mode of 
strings connecting these branes becomes massive since the mass square of this string 
is given as L 2 / R A — l/(27ra') [31]. However, we must notice near the mid point of the 
baryonium configuration where the two branes are overlapping, then we might expect 
to find tachyons at this point. Here, we study this problem through the fluctuation 
mode around whole region of the classical configurations. 





(17) 
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3.1 Fluctuations 



In order to consider the fluctuations, we back to the action of D5 brane (jlj), and 
expand it with respect to the fluctuations, Sr(t, 0) = r — f, 8x(t, 0) = x — x and 
8A t (9,t) = A t — A t , up to their quadratic terms. Here f, A t (9) and x are the solutions 
of the equations of motion. As for the SA t , we retain it although it is not dynamical 
due to the lack of the kinetic term. The quadratic part of (pO) for these fluctuations is 
given as 

L (a) = -sin*0 (L% + Lfa) , (18) 



where 



B, 



(0) 



8r' 2 + {^) 5x' 2 - (r 2 + r' 2 ) 5x 2 - + x'A 8r 2 + 2x'r'8x5f 



+d r (B {0) ) r'5r'5r + d r [ ^4#(o) I x'5x'5r 



2G (0 ) V Vfl/ 



r 5r + „ x + 2 „ r'x'Sr'Sx' , (19) 



(20) 

where r' = dor {6), x' = dox{6), 5r = d t 8r and x = d t 8x. 

G ( o)=(r 2 + (0 2 + ^(^) 2 )^. (21) 

B {0) = e**/y/§^ , A {2) = \d 2 r (e*/ 2 y/S^j , e* = 1 + ± . (22) 

We notice here as mentioned above that 5A t is not a dynamical fluctuation since 
there is no time derivative term for it. It can be regarded as a kind of an auxiliary 
field, then we integrate it out from L( 2 ). It is performed by rewriting Lif£ as follows, 

L* = -C (2) (5A> - yjl* ind /C (2) ) 2 + L* ™ , C m = (0) , (23) 
where 

(24) 
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In order to remove 5A t , we replace L,^ by L^t md in (|T8j) . then we can restrict to two 
fluctuations 5x and 5r to see the frequency of these fluctuations, which have their time 
derivatives. 

The modified quadratic term is obtained as 



j — jx,r , j At ind 
M2) = L {2) + L {2) 



(25) 



and 



\2) 



+ (r 2 + r' 2 ) (-Si 2 4 I - 



K , /2 
~V + X 



-5r 2 + {—) -J— 5r' 2 ) + 2x'r'5x5r 



R) Q(o) 



R) Q(o) 
+Q(i)r'5r'5r + Q^ 2 )x'5x'Sr 



-5x' 2 



-2x'r'Sx'5r' 



RJ Q(o) 



where 



A(2) = A(2) + 



1 e 2 *G (0) D 2 

2 -B(o) Vv 

D 



(0) 



2 , Q(0) = r 2 + r' 2 + x' 2 (^) 4 



(26) 



(27) 



Q(i) = ( <9 r (5(0)) - ^e*9 r (e *B (0) ) J , Q( 2 ) = ( d r 



V 4 5, 



(0) 



r A D 2 



R 4 



(28) 

Then we can estimate the eigenvalues of the frequency of fluctuation by imposing an 
appropriate boundary conditions for each eigenfunctions for these fluctuations. These 
eigen functions are obtained solution of 



( M„ 


M rx \ 




1 M xr 


M xx j 





where 

M„ 

M xx 
M rx 



1 

2 [r 2 
+Q(i)r'd e 

j~2 _|_ ^,/2 



where 5/j = (8r, 5x) 



V 4 sin 4 9 B {0) (R A 



(29) 



+ *' 2 ) a (co a? + i (2) - -^d 9 f ( ^ + x' 2 1 ^ 



sm^ "V 2i? 4 Q (0 ) V^ 2 



1 /r 4 sin 4 O (0) 

On 1 



1 \ Q(2)x'd e + -r\^d t 

2 [ sin 



sin 4 9 V 2i?4 Q(o) 
r 4 sin 4 6 



r 2 + r' 2 



\R 4 Q 



(0) 



(30) 
(31) 

(32) 



1 



1 a ■ 4 an i , 1 a ^ r4 sin 4 6» , , 



sin 4 9 V R4 Q 



(0) 



2\ sin 4 fl' 

-rVB (0) ^} , (33) 

In the followings, we estimate the eigen-frequencies of the fluctuations in order to see 
whether the solution is stable or not. 
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3.2 Estimation of eigen-frequencies 

By assuming the following form for the fluctuations, 

5r(t, 6) = e iuirt <f) r {6), 5x{t, 6) = e**-*^) , (34) 

we estimate the value of frequency u rjX by solving the equations given in (1291) for 
$o < & < 7T- Since we are intersted in the lowest eigenvalues of u; ria; to see the stability 
of our solution. Such fluctuations are belonging to the no-node wavefunctions for 
4> r ,x{9)- Then the equations are solved under the condition, (fi r ,x(@o) — const, and 

<M?r) = (35) 

or 

de<fi r ,x( n ) = (36) 
for Diriclet or Neumann condition at the boundary 9 = ir respectively. The constant 
values at 8q are determined by the normalization condition of the wave functions 4> r ,x{d)- 
We calculated the eigenvalues for the classical configurations, which provide the 
minimum value of U for each u, since such configurations would be stable and positive 
u 2 s are expected. In the Fig. [31 we show the lowest eigenvalues for the Diriclet 

1.0 r 

0.8 - 
0.6 - 

0.4- . 
0.2 - 

• 

•••• • . , * 

00 [ — , — , — , — i — , — , — , — i — , — , — , — i — , — , — , — i — , — , — , — i v 
0.0 0.2 0.4 0.6 0.8 1.0 

Fig. 3: The lowest for Diriclet boundary condition. The upper (lower) points show 

condition for various v. For all cases, we can see the baryonium configurations are 
stable against the fluctuations. 

As for the Neumann condition, the lowest mode of 5x is the zero mode (u x = 0) 
with 4> x (9) = const., so we show the lowest eigenvalues of u 2 in the Fig. HJ Again, we 
find positive u 2 for any v which implies the stability of our solutions for any v. 

We notice here the following points. 
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Fig. 4: The lowest u 2 for Neumann boundary condition. 



(i) Firstly, the lowest eigenvalues are obtained for the zero-node wave- functions of 
each fluctuations, and the values of u 2 become large with increasing node number of 
the wave-functions as expected. They are abbreviated here. 

(ii) In the second, we show above the eigenvalues at U — U min , however, we could 
find positive eigenvalues at other values of U. In this sense, all the classical solutions 
obtained here are stable against small fluctuations around those configurations. 

(iii) The next point to be noticed is that the above analysis is given for < v < 1, 
and the point v — and 1 should be excluded since the configurations in these limit 
of v are not the baryonium state. 

3.3 Expected configuration for N c = 3 

The baryoniums are composed of the vertex, n quarks and n anti-quarks. This combi- 
nation of quarks is denoted as (n, n), As a practical matter, we consider here the case of 
N c = 3, then we find three baryonium configurations of (3, 3), (2, 2) and (1, 1) quarks. 
The one of (2, 2) and (1,1) are shown in (a) and in (b) of the Fig. [5] respectively. They 
will give a good guide for baryonium hunting in the future. 




5- X 

(a) Cb) 



Fig. 5: Expected two baryonium configurations for N c = 3 are shown. 
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4 Stability of D8/D8 and D6/D6 as flavor branes 



In the case of type IIA string theory, a holographic gauge theory is considered in the 
D4 brane stacked background, 

//7\ 3 / 2 /ff\ 3 / 2 / dJJ 2 \ 

ds 2 10 = ( 1 ) (v»vdx»dx» + f(U)dr 2 ) + (-J + U 2 dQlj , (37) 

where 

m = i _ (^) 3 , (38) 

and Ukk is related to the compact size of S 1 in the D4 world volume. Then the flavor 
degrees of freedom has been introduced by embedding D8/D8 [30] and D6/D6 [28] 
probe branes. 

In these cases, the fundamental strings are not dissolved in the bound state branes. 
So the fluctuations are not affected by them as in the baryonium configurations. Our 
purpose to study these cases is to see the stability of those configurations agaist the 
fluctuations by using our method given above. 

While, for a special case in [30J, the stability of the configuration has been shown as 
mentioned below, we can show their stability for more general configurations by our 
method. In [28] , the proof of the stability of D6 / DQ is abbreviated, but its stability is 
assured in the followings. 



4.1 DS/D8 branes 

For the case of D8/D8 probe, they are embedded with the following induced metric, 



(39) 

The embedded configuration of the probe brane is obtained from the D8 brane action, 



^3/4 



S D8 oc J d 4 xdUU\U/R) 
where r' = dr/dlf. The equation of motion for r(U) is obtained as 



V (c,)(t ' )2 + (f) 3 7M' (40) 



t' 



U\U/RW* . = co , (41) 

^) 2 + (#) 7m 

where Co denotes an integral constant, which is taken here as 

c Q = U^f 1 / 2 (U )(U /R)^. (42) 
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Here Uq(> Ukk) denotes the minimum value of U for the solution. Then the solution 
is given as 

r u /R\ 3/2 1 1 
t(U) = ±c dU(£) ..... . (43) 

ho \UJ f(U) ^/u 8 f(U) - cl 

As in the case of the baryonium (see Eq. (fl6l) ). also this solution can be interpreted as 
the bound state of D8 and anti-D8 branes which are connected at U = Uq. Then we 
are able to study the stability of this configuration as above by checking the existence 
of the tachyon on the brane. The procedure is parallel to the case of the baryonium. 
By setting the fluctuation as St = r — f , where f denotes the classical solution given 
by Eq.(j4"3~|). the D8 action is expanded by St as, 



{R/Uf 



2GIP(U) 

rTd^rdjT + ■••), (44) 



2G T f(U) 
where 

G T = r' 2 + (R/Uf/f 2 (U) . (45) 

While, at this stage, we can see that there is no tachyon for any Uq from this form, we 
proceed the above procedure. By setting as St = f T (x^)(f> T (U), the eigenmass equation 
of St is obtained as follows 



where ml is defined as 



T 



rrWvUi??) = m 2 J T (x») . (47) 



As in the previous section, we can investigate numerically the eigen mass m 2 . of the 
fluctuation <p T thorugh (|46l) . In this case, however, the boundary condition of <p T is 
taken as 

hm <f> T (U)cxU-i, 7 >^ (48) 

in order to have a normalizable wave function. This is assured from the Eq. (I46p . Under 
this condition, we find positive value as the lowest mass eigenvalue m 2 for any Uq. The 
results are shown in the Fig. |6l So we could show the stability of this brane anti-brane 
bound-state also. We should notice that m T increases linearly with Uq/Ukk as seen 
from Figj6j 

Finally we notice that the stability of this configuration has been studied for the 
case Uq = Ukk [30J . However, here, we could show its stability for any value of U , 
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25 r 




n=Q 



Fig. 6: The lowest value of m T versus Uq/Ukk for D8-D8. The values are for zero node 
wave- functions, n = 0. 



Uq > Ukk- The reason why this bound state is stable for any Uq, where D8 and anti- 
D8 branes are connected, is that there is no U(l) flux in the brane in this case. This 
phenomenon seems to be universal for any bound state of D and anti-D branes. 



4.2 D4/D6-DQ branes 

Instead of ([3"T]) . the D4 background solution is given in [28] as follows 

/ TT\ 3/2 

ds\ Q = I - 1 (p„ u dx^dx v + f{U)dr 2 ) + K{p) (dX 2 + X 2 dQ 2 + dr 2 + r 2 d<j) 2 ) . (49) 



U 



3 \ 2/3 



where p 2 = X 2 + r 2 and 

The D6 is embedded in the following metric, 



R^U 1 / 2 



(50) 



(51) 



/7T\ 3/2 

ds 2 = (-) ri^dx^dx" + R 3/2 U 1/2 dtt 2 2 + 



RJ 



(i) 3/2 w) 



dU 2 . 



(52) 

Here we supposed that the profile is determined by t(U). Then the DBI action of D6 
brane is given by 

^17/4 



W6 — J-D& 



T m R 3/2 



R 
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f(U)G T . 



(53) 



where G T is given in Eq.( l45p above. Solving the equation of motion derived from L^q, 
we obtain DA/ D6-D6 embedded solution t(U) in the background ( H9l) as in the D8 
case. The equation and its solution are given as follows, 

(U/RY^f^^= = d , (54) 

where the integral constant do is given as 

d = f 1 /\U )(U /R) 17 / 4 . (55) 
Here Uq is the minimum value of U and t'(Uq) = oo, and the solution is given as 




Fig. 7: The lowest value of m T versus Uq/Ukk for D6-D6. The values are for zero node 
wave- functions, n = 0. 



As for the fluctuation of r, we can see its quadratic term by expanding (|53|) with 
respect to 5r(— r — r c i) as follows, 



S m oc J d 4 xdU(U/RY 7 /^f(U)G T \l + 2 ffi^ m (Sr') 2 + 
{R/Uf 



2Glf\U) 

2G T f(u/ Vd ^ 6r + --}> < 57 > 
By setting as 5r = f T (x^){p T (U), the eigenmass equation of <p T (U) is obtained as follows 

5 /4 \ . (TT/R\5/4 } 



-{U/Rfl\ \ / ,, rn „ (U/R) 



du l y /f(U)G T ^j^ d u j + ^f(U)G T ^ nu) cp T (U) = , (58) 
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where the mass eigen- value m 2 is defined as above, f]^ y d^d u f T {x^) = m^frlx 11 ) . The 
linear term of 5t disappears due to the field equation of r. 

In solving ( l58l) . we must impose the boundary condition for if from the normaliz- 
ability as 

lim <p T (U) oc If-"* 1 , 7! > I . (59) 

U^oo 8 

Under this condition, we have find positive mass eigen-value m 2 as in the case of D8 
branes. The results are shown in the FigJTl In this case also, m T increases linearly with 
Uq/Ukk as seen in the case of D8 branes. In any case, the D6/D6 state is also stable 
against the fluctuation of r(U). 



5 Summary 

We have examined the stability of D brane anti-D brane bound state embedded as 
a probe in some 10D supergravity background. We have restricted our attention to 
some special configurations, the baryonium state in type IIB theory and flavor branes 
in type IIA theory. 

The vertex of the baryonium has been obtained as a bound state of D5 and anti-D5 
brane in type IIB theory in terms of the D5 brane action, which wrapped on S 5 in the 
10D AdSs x S 5 and gives baryon vertex at the same time. The configuration of the 
baryonium is determined by the main two parameters, the gauge condensate q, which 
would determine the enrgy scale, and the flux number (1 — v)N c (or vN c ), which is 
proportional to the quark number in the baryonium. 

Previously, we could find an appropriate solutions of the baryonium which gives a 
minimum of the classical action for any value of v and q. In order to confirm the stability 
of this configuration furthermore, this time, we have studied the time-frequency u of the 
fluctuations for the relevant fields living on the D5 brane to search for the unexpected 
tachyonic modes. The eigen value cu has been investigated numerically for the eigen 
functions (j) T {9) and (j) x {d) on which the Dirichlet or Neumann boundary conditions 
are imposed. For both bondary conditions and at any value of u, we have found non- 
negative oj 2 for the configurations given at U = U min . Surprisingly, the non-negative uj 2 
is observed even if we consider the fluctuations around any configurations of U(> U m i n ). 
In this sense, any baryonium configuration given as the classical solution would be 
stable agaist the fluctuations around the baryonium solutions considered here. 

As in the case of the baryon, we must add quarks to the vertex to find a complete 
baryonium state. The method would be parallel to the case of the baryon [18], where 
the no force conditions are needed at each cusp to connect (1 — v)N c fundamental 
strings, this will be done in the near future. 

We have also investigated the stability of the bound-state of D8/D8 and D6/DQ, 
which are introduced as the brane for the fundamental quarks in the type IIA theory. 
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We found that these bound states are stable. This stability would be promised by the 
fact that there are no U(l) flux dissolved in the branes in these cases. 
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Appendix 

A: Comment on Equation of motion 

And the action is rewritten by eliminating the gauge field in terms of the solution for F t g 
given by the second equation of fflOl) to obtain an energy functional of the embedding 
coordinate only0: 

U = -4-; / d0 e* /2 Jr 2 +r /2 + (r/RYx 12 JVJO) . (60) 

V u (9) = D{u,9) 2 + sin 8 9 (61) 
where we used T 5 fi 4 i? 4 = N/(3ii 2 a r ). 

We notice that the solution for x(9) in our baryonium configuration is given by a 
two valued functions for the variable 9, which is restricted to (i) 9 < 9 < n or (ii) 
< 9 < 9i as explained below. So it is convenient to solve the equations by changing 
the variable from 9 to x in U as, 

U = -^7 / dx e * /2 ^r 2 9 2 + r 2 + {r/Ry ^/vj9), (62) 

where dots denotes the derivative with respect to x. In this form, x is not contained 
explicitly, then we can introduce an integral constant h as a "Hamiltonian" for the 
corresponding "time variable" x as follows 

h — rp r + 9pg — L , (63) 



where L = e*l 2 Jr 2 9 2 + r 2 + (r/i?) 4 JV U (9) , 



dL dL 2A„ „ f R \ 2 

Pr= dr~ = rQ ' Pd = ~d9 = = Q ' Q= W\ 



1 U is obtained by a Legendre transformation of L, which is defined as S = J dtL, as U 
Tr£-Fto — L. Then equations of motion of © provides the same solutions of the one of U. 

dF te 
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Then h is rewritten in terms of the momentum as 

h = ■ 



RJ \ 



(65) 



and the Hamilton equations of motion are obtained as 



r= Q' 6 



dh 



dh 



r 2 Q ' ^ r dr ' ^ 9 86 



(66) 



These equations are convenient to find the baryonium vertex solutions. 
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